りこ の ま キ 


(し Ez 


工 


SSN 


ここ で は , テー ブル を 用 いて 除算 回 路 を 実現 する 方 法 を 解説 す 
る . テー ブル を 使用 し て も , 補正 を 加え る こと で ある 程度 の 精 
度 を 維持 する こと が で きる . な お , 本 稿 で 設計 し た 回 路 の ソー 
ス ・ コ ー ド は 本 誌 付属 CD-ROM に 収録 し て いる . (編集 部 ) 


除算 回 路 を 過去 4 回 に わた っ て 説明 し て きま し た が , 今 
回 で 本 連載 も 終わ り で す . 
連載 第 1 回 本 誌 2005 年 8 月 号 , pp.109117) と 第 2 回 同 
年 10 月 号 , pp.136-144) で 述べ た 回 復 法 と 非 回 復 法 は , け 
た 単位 で 商 を 確定 し て いく 方 法 で し た . また , 第 3 回 同年 
11 月 号 , pp.100-109) で は その 高速 化 に つい て 述べ まし た . 
第 4 回 同年 12 月 号 , pp132-138) は 漂 近 法 に より , けた を 
超え て 真 の 商 に 近づけ る 方 法 を 紹介 し まし た . いずれ ぞ 繰 
り 返し 演算 が 必要 で ある 」 と いう 点 で 共通 し て いま し た . 

今回 は テー ブル 近似 に より ,「 繰り 返し 演算 を 伴わ ず に 
商 に 直接 た どり 着く 」 方 法 を 取り 上 げ ま し ょ う . テー ブル = テ 
メモ リ と いう 割り 当て を 考え た と き , メモ リ が 潤沢 で 使い 
や すく な っ た 近年 の FPGA に お いて 非常 に 有効 な 方法 で あ 
る と と も に , 除算 だ け で は な く 初等 超越 関数 三角 関数 や 
対数 関数 ) の 実現 に も 応用 で きる の で , ぜひ と も 活用 し て 
いた だ きた いと 思い ます . 

「 テー ブル で は 精度 が 今ひとつ 」 と いう 批判 を 聞い た こと 
が あり ます . お そら く 商 を すべ て テー ブル 索引 だ け で 処理 
し よう と する か ら で し ょ う が , それ で は あま り に 芸 が あり 
ませ ん . 各種 補正 を 加え る と いっ た く ふ う を 施す こと で , 
精度 を 上 げ ら れ ま す . 


人 @ 除算 で は な く 逆 数 関数 で 実現 
あえ て 述べ る まで も な いこ と で す が , 除算 は 除 数 と 被 除 
数 と いう 二 つ の 引き 数 が 必要 と な る 演算 で す . テー ブル を 
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テー ブル 近似 と 精度 向上 の た め の 補正 


関連 デー タ 


- 鈴木 昌治 


EE 、 


用 いて 除算 を 実現 する 場合 で , 二 つ の 引き 数 の 組み 合わ せ 
と な る と , どう し て も テー ブル が 肥大 化し て し まい , 現実 
的 で は あり ませ ん . この あたり の 事情 は , 前 回 述べ た New- 
ton-Raphson 法 と 似 て いま す . そこ で , 今回 も Newton- 
Raphson 法 の 場合 と 同じ く , 逆数 関数 を 実現 する こと に し 
ます 、 

連載 第 2 回 の コラ ム で ,「 除算 は 浮動 小数 点 と 相性 が 良い 
演算 で ある 」 と 述べ まし た が , ここ で は その 特性 を 利用 し 
まし ょ う . と いう の も , か り に 入力 で ある 分 母 7 を 正 数 に 
限っ た と し て も , すべ て の 範囲 を テー ブル で カバ ー す る の 
は 不可 能 だ か ら で す . 有限 な テー ブル 上 で 効率 良く 関数 を 
表現 する に は , 逆数 の 場合 , 浮動 小数 点 表現 を 用 いて 次 式 
の 範囲 に 入力 を 限定 する こと で す . 


これ に より 分解 能 を 決定 すれ ば , テー ブル の 深 さ ( アド 
レス 幅 ) を 決め られ ます . また , 式 1) に よっ て 商 ? は , 式 
( 2) の よう に な り ま す . 


要求 され る ビッ ト 精度 が 決ま れ ば , テー ブル の デー タ 幅 
を 決定 で きま す . 

以上 で , テー ブル を 用 いた 逆数 演算 の 土台 が で き 上 が り 
まし だ 、。 


" 基本 は 直線 補間 


の 9 と も に 8 ビット 程度 の 精度 が 要求 され る 場合 , 2K ビ 
ッ ト の メモ リ が 必要 と な り ま す . 今どき は , すべ て テー ブ 
ル 化 し て し まっ た ほう が 現実 的 で し ょ う . し か し , 16 ビ ッ 


5 ュ で トー /。 リュ 
@ C に 6 C て に = | [) 宮 口 | 2 AU ー 
一 ーー ーー 一 ーーー ドー ーー ーー ニー 2 
C C 6 C に @ C C ぐ C @ Cc ご と て て 4 
(し ーーーー 


2 5 の 計 ご の 置 * 

ト の 精度 が 入出 力 に 求め られ る 場合 に は , メモ リ は 1M ビ { 区 
ッ ト と な り , ASIC や FPGA に 内 蔵 す る テー ブル と し て は 関数 値 o) 図 
大 き な も の に な り ま や 除算 の 前 半 部 で し か な い 逆 数 演算 10 

に 許さ れる サイ ズ で は な い ). 還 症 

ここ で 登場 する の が 直線 補間 に よる 補正 方 法 で す . 定義 補 Es し 
域 を いく つか の 区 間 に 分 割 し , その 境界 に お ける 関数 値 を 《 1+ 上 
テー ブル 上 に 持ち , 一 つの 区 間 の 値 を 直線 補間 に よる 演算 了 
で 求め る 方法 で す . 

7 

@ 直線 補間 を 用 いる と 回 路 は シン プル , 精度 は いま いち 較 間 

図 1 に , 直線 補間 の 原理 を 示し ます . ここ で は 見 や すく 1 1.5 20 
する た め に 二 つ の 区 間 に 分 割 し , 1 ミ / ぐ 15 の 区 間 を 直線 図 1 直線 補間 に よる 逆数 演算 の 原理 


1 ミ 9 ぐ 2 の 正規 化 区 間 を さら に 分 割 し , 各区 間 に 直 線 補間 を 適用 する . 


補間 し て いる よう す を 表し て いま す . 求め る 値 を 2=1 十 ん 
の と き の 値 と すれ ば , 7( 1 十 ん ) は 式 3) で 与え られ ます . 


コキ 人) ユー の ※ メ た リオ た 15) の ot ( 3) に た と 還 ECHOOPREPEEERREEEEEEEEREEEREERERE 《 4) 
まず は , 入出 力 16 ビ ッ ト と し て 話 を 進め まし ょ う . 図 1 出力 の 16 ビ ッ ト 精度 を 目標 と し て いる の で , テー ブル 上 
の よう に 葉 1) の 区 間 を 2 分 割 と いう わけ に は いか な い の の 境界 点 デ ー タ は 2 ビッ ト ほど 余裕 を みて 18 ビ ッ ト と し ま 
で , ここ で は 232 分 割 す る こと に し ます . テー ブル の アド す . し た が っ て , 32 区 間 X 18 ビ ッ ト xX% 区 間 の 始点 と 終 
レス は 5 ビッ ト に な り ま す . 式 ) よ り 42 の MSR most 点 の 2 デー タ が 必要 ) で 1K ビ ッ ト 強 の テー ブル と な り ま す . 
significant bit) は つね に 1 な の で , これ を 除外 し た 小 テー ブル か ら 出力 され た 区 間 の 始点 と 終点 デー タ を , 
数 点 以下 の 部 分 を 16 ビ ッ ト に 割り 当て ます . その た め , d 10: 0{ テー ブル を 引く の に 使用 し な か っ た 部 分 ) で 直 
[ 15: 11] が テー ブル を 引く た め の 5 ピ ッ ト に な り ま す . こ 線 補間 し , 丸め を 行え ば 演算 完了 で す . 
こ で , 1 区 間 の 幅 を ヵ と する と 次 式 の よう に な り ま す . ブロ ッ ク 図 を 図 2 に , 演算 の けた の 割り 当て を 図 3 に 示 


関数 実現 の た め の 直 線 補間 


本 稿 の よう な 初等 超越 関数 を 実現 する 際 に 直線 補間 を 用 いる 場合 , 関数 な の で 差分 は 減算 され る が , 増加 する 関数 で あれ ば 加算 と な る 
本 文 の 図 2 に 示し た 構成 を と る 必要 は あり ませ ん 本 文中 で は , あく 差分 を 符号 付き と すれ ば 加算 で 統一 し て も か まわ な い ). 
まで 原理 的 表現 と し て ある ). その 理由 は , 直線 補間 が , 方 , 後者 の 理由 か ら ぱ 前 者 と も 関連 が ある が ), 差分 テー ブル 
e 補間 する 2 点 が テー ブル で 与え られ る は 終点 テー ブル に 比べ て デー タ 幅 を 小さ く で きま す . 本 文 で 示し た 
e 定義 域 が 限定 され た 場合 , 変化 率 が 穏やか で ある 例 で は , 終点 テー ブル は 18 ビ ッ ト 幅 で し た が , 差分 テー ブル で あれ 
と いう 条件 を 持つ か ら で す . ば 12 ビッ ト 幅 で すむ た め , テー ブル ・ サ イズ を 小さ く で きま す . 当 
前 者 の 理由 か ら , 無 作為 な 同じ ディ メン ジョ ン を 持つ 2 変数 問 を 補 乗算 器 も 小さ い 規 模 の も の で すみ ます . 
間 す る の と 異な り , 直線 補間 で は その 区 間 の 始点 と 次 の 区 間 と の 差 これ に よっ て , 図 2 は 図 A-1 の よう に 表現 で きま す . 
分 を 与え られ れ ば 補間 可能 と いう こと に な り ま ず 始点 や 終点 は 位置 , 


差分 は 距離 と いう ディ メン ジョ ン の 違い が ある ). テー ブル の 中 身 は 始点 図 

何で も 可能 で すし , 関数 その も の は 既知 の も の な の で , か な ら ず し 15: Sa | 一 g 
も 始点 と 終点 を 与え る 必要 は あり ませ ん . これ に より , 始点 を g, 差 
分 関数 を # *) と する と , 本 文中 の 式 3) は 次 の よう に 表現 で きま す . 


PA g 十 9 た ヨコ 4 Z)ー た X & Z) チキ で サザ ド すす さす キヤ で そ ヤ テ ザ ャ ャ ャ マヤ ( A-1) 


明らか に 演算 が 簡略 化 さ れ て いま ず 逆数 の 場合 は 一 様 に 減少 する 図 A-1 直線 補間 の 別 の 表現 方 法 
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十 X ち 


区 間 終 点 図 終点 図 
テー ブル 図 


2 直線 補間 に よる 逆数 演算 ブロ ッ ク 


9 の 上 位 け た で テー ブル を 引き , 得 ら れ た 区 間 の 始点 と 終点 を 9 の 下位 の 
けた で 補間 する . 


k 


リス ト 1 図 2 の 回 路 の RTL 記 述 


modu1e div tab1e a ( 
す  , // quotient [16] :inE [15:00] :Eraction 
d ): // divisor [15:00] : Eraotion ( [16] : static nt 1) 


// 1o 
1npu も [15:00] d : 
outpu セ も [16:00] q : 


// Eunction (演算 処理 部 ) 
// tab1e 
Funot1on [35:00] Etab1e : 
// テー ブル [35 :18] : 始 点 [17:00] : 終 点 
ュ 1npu も [15:11] d 2 
で a8@e 
Ftab1e 18'h20000, !h1F07C 
Ftab1e 18'h1EF07C, !h1E1E2 
Ftab1e 18 !h1E1E2 , !h1D41D 


Ftab1e !h10C97, !h10842 
Ftab1e !h10842, !h10410 
Ftab1e !h10410, !h10000 


endoa8e 
endFuno も tion 


Wire [17:00]  p stt 
Wire [17:00] p end 


as81gn {p_stt,p_end) = ftab1e (d[15: 


// bi-1inear interpolatron (直線 補間 ) 
wire [29:00] inErp = ( 12'h800 - d[10:00] ) * p st 上 


+ dQ[10:00] * p end + 12'h800 : 


a881gn q = ntrp[28:12] : 


endmodu1e 


し ます . また , この 回 路 の RTL 記 述 を リス ト 1 に 示し ま 
す . 出力 の MSB は 2= 1 の と き 以 外 は つね に 0 な の で , 
実質 の 16 ビ ッ ト と いう 意味 で 出力 は 17 ビ ッ ト と し て あり 
まず 原理 に 忠実 に 表現 し て いる の で , 実際 の 構成 に つい 
て は , p.129 の コラ ノル 関数 実現 の た め の 直 線 補間 」 を 参照 ). 

2 を 見 て わか る の は , 非常 に シン プル な 構成 で ある と 
いう こと で す . これ が , 直線 補間 を 用 いる 最大 の メリ ッ ト 
で し ょ う . 反面 , 精度 は それ ほど 高く な いと いう デメ リッ 
ト も あり ます . 誤差 を 薄 5) の よう に 定義 し ます . 


eg テル 9 の )- を の) 8 1108ieie25 ( 5) 


定義 域 に お ける 誤差 の よう す を 図 4 に 示し ます . 各区 間 


区 
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16 12 8 4 O 
gh 
16 12 8 4 0 
{g), (ag+ 応 図 
16: 12 8 4 0 
9 
図 3 直線 補間 に よる 逆数 演算 の けた 割り 当て 


を る と を a 十 尺 は , それ ぞ れ テー ブル 出力 の 始点 と 終点 の デー タ を 表す . 


0.00025 


o00020 | 


0.00015 


誤差 較 


0.00010 


0.00005 


0 
1.0 1.2 1.4 1.6 1.8 2.0 


図 4 直線 補間 に よる 逆数 演算 の 誤差 
9 が 1 に 近い ほど 誤差 が 大 きい の は , より 曲率 が 大 きい ( 直線 に な じみ 
に くい ) か ら で あ る . 


の 中 央 に 誤差 の ピー ク が あり , それ が 分 割 し た 区 間 の 数 だ 
け 並 ん で いま す . 2 が 1 に 近い ほど 誤差 が 大 きい の は , よ 
り 曲率 が 大 きい 直線 に な じみ に くい ) か ら で す . 

結果 的 に は , 16 ビ ッ ト 精度 を 目指 し な が ら も 12 ビ ッ ト 
程度 の 精度 し か 得 ら れ て いな いこ と が わか り ま す . ここ で 
注意 し て ほし い の は , 図 4 は RTL の シミ ュ レ ーション 結果 
を 表し て いる わけ で は な く , 米国 Microsoft 社 の Excel で 
シミ ュ レ ーション し た も の だ と いう こと で す . Excel で は 
10 進 で 15 け た の 精度 が ある と いう こと な の で , テー ブル 
や 補間 演算 な ど は か な り の 精度 で 演算 され て いま す . つま 
り , 直線 補間 を 用 いる 限り , テー ブル や 演算 の 精度 を 上 げ 
た と ころ で , これ 以上 は 望め な いわ け で す . で は , どう す 
れ ば 精度 を 上 げ ら れる の で し ょ うか . 


人 @ 区 間 を 細分 化し て 精度 を 上 げ る 

1 か ら ,「 区 間 を も っ と 細か く すれ ば 精度 を 上 げ ら れ 
る 」 と いう こと が 直感 的 に 理解 で きる と 思い ます . し か し , 
や みく も に 試行 錯誤 で 精度 を 上 げ て いく の は , スマ ー ト な 
方 法 と は いえ ませ ん . まず , 区 間 の 分 割 数 と 精度 の 関係 を 


ーーー ーー ーーー 一 | 王 デ トド ーー ーー 
て い し 1 に 介 に に 6 介 
ーー 一 | 一 ーー ーー トー し ーーー 
と に て @ て 


図 5 

区 間 分 割 の 違い に 
よる 誤差 の 推定 

区 間 を 細分 化す る こ 


と で 精度 を 向上 で き 、 ! 
る . 4 ag 填 万 4 ag 填 /2 g 十 


考察 し て み ま し ょ う . 

図 5 に 区 間 分 割 の 違い に よる 誤差 の 違い を 示し ます . 区 
間 の 始点 を 区 間 の 幅 を ヵ と し ます . 誤差 の 最大 値 が 区 
間 の 中 央 に ある と 仮定 する と き , 最大 誤差 は, 式 6) で 
与え られ ます . 


8i 全 5 ( 6) 


1 1 1 
ーー 2 
昭 g 十 刀 /2 ーー ル 


区 間 の 幅 が ヵ /2 と な っ た 場合 の 最大 誤差 は , 再 7) で 与 
えら れ ま す . 


1 1 1 
ー 十 / 2 の に 2887 7 
o+ が 7 4 : ーッ / 


以上 より , 区 間 が 1/2 と な っ た 場 洪 区 間 の 分 割 数 を 2 倍 
と し た 場合 ) の 誤差 の 比 は , 次 式 で 表 さ れ ま す . 


誤差 が 最大 と な る 区 間 は 7= テ 1 を 始点 と する 最初 の 区 間 
で ある こと は わか っ て いま す . そこ で 。ー1 と し て , ヵ が 分 
割 を 進め る ほど 微小 区 間 と な る こと を 考慮 し て , ヵ ー0 の 
極限 を と り ま す . 


。 2 」. カカ 二 1 1 
hm = hm ニ 
っ 0g] 4 っ 0 カカ +4 4 


hie763918 の SS289 Sie ( 9) 


葉 9) か ら , 区 間 の 分 割 数 を 2 倍 に する と 精度 は 4 倍 に な 
る こと が わか り ま す . 実際 に は ヵ = テ 0 で は な い の で , 4 倍 よ 
り は 若干 小さ な 値 に な がり ま す が , お よそ の 目安 に は な る で 
し まう 。 


に 
ーー 0 ーー 一 
“ ンー ビ e 遇 7 みーこ 
Go ゎ (OO ーー 2 ・ で = 
0.00025 | 
0.00020 
図 0.00015 
押 
器 
0.00010 
0.00005 
0 
1.00 1.01 1.02 1.03 
び 


図 6 直線 補間 に よる 最初 の 区 間 の 誤差 
誤差 の ピー ク が 中 央 に あっ て , ほぼ 左右 均等 の 誤差 た 示し て いる こと か ら , 
偶数 次 成分 が 大 きい こと が わか る . 


正規 化 2 次 誤差 関数 較 


補間 直線 か ら 減算 較 
する こと で , 真 値 に 
近づけ る 図 


と ここ 


- 
正規 化 2 次 誤差 関数 芝 
y4% 1- x) 図 
図 7 ( 
2 次 曲線 に よる 誤 1「---- ドーー ン | 
差 補正 
真 値 と 直線 補間 の 差 


( 誤差 ) を 2 次 関数 で 
補正 する . 区 間 始 点 =0 補間 位置 較 


区 間 終 点 三 1 


先ほど の 32 分 割 で は 12 ビ ッ ト の 精度 だ っ た の で , 16 ビ 
ッ ト に 追い 込む に は 4 倍 の 128 分 割 と すれ ば よい こと が わ 
か り ま す . この 程度 で あれ ば 許さ れる テー ブル ・ サ イズ で 
す が , さら に 20 ビッ ト , 4 ビッ ト と 高い 精度 が 必要 と さ 
れ た 場合 に は 規模 的 に 厳し く な っ て きま す . また , 精度 を 
上 げ る の で , デー タ 幅 も 大 きく と ら な く て は な ら な いこ と 
が , これ に 拍車 を か ける こと に な り ま す . 分 割 数 を 上げる 
だ け で は お の ず と 限界 が 見 えて きま す . 


[2 高 次 補正 を 採用 し て 精度 を 上 げ る 


直線 補間 を 級数 展開 な ど に お ける 1 次 まで の 近似 と と ら 
える な ら ば, 2 次 的 な 補正 を 加え る こと で 精度 を 上げ ら れ 
る と いう こと が 考え られ ます . 初等 超越 関数 を 級数 展開 し 
た 場合 , 高 次 に な る ほど 係数 が 小さ く な る の が 一 般 的 な の 


Design Wave Magozine 2006 Jonuary 131 


で , 「 1 次 の 次 は 2 次 」 と いう の は きわ め て 
あの 生き 


自然 な 考え か た 


信 2 次 補正 で 精度 を 20 ビッ ト に 向上 

こ 直 線 補間 32 分 割 ) の 場合 の 最初 の 区 間 に お ける 誤 
差 を 拡大 し て 示し ます . 誤差 の ピー ク が 中 央 に あっ て , ほ 
ほぼ 左右 均等 の 誤差 を 示し て いる こと か ら , 偶数 次 成分 が 大 
きい こと が わか り ま す . この 形 で あれ ば , 2 次 の 補正 を か 
けれ ば 精度 を 向上 させ られ る こと が 予想 で きま す . そこ で , 


款 叶 に 半 当 が tee0eNRRAOWHKeeeiBdtee teeeesnietet 《 10) 


と いう 2 次 関数 を 定義 し , 分 割 さ れ た 1 区 間 を 0 ミ 7 ぐ 1 に 
割り 当て まず 図 7). つま り , 直線 補間 に 用 いた d 10: 0O] 


れ た 2 次 誤差 補正 関数 と 考え , 
直線 祖 


ーション 結果 ). 


ト 2 に 示し ます . 


この 


図 11 に , 最初 の 区 間 に 
図 か ら わ か る こと は , 1 区 間 内 の 誤 
よっ て 占め られ て いる こと で す . 


区 間 ご と 


お ける 誤差 を 


れ た こと に 


に ゲイ ン を 変え て 


和 間 値 か ら 減算 し て や れ ば よい の で す . 2 次 誤差 補正 
を 施し た 後 の 誤差 を 図 8 に 示し まず Excel( 
この 時 点 で の 精度 は 20 ビ ッ ト 静 ま で 
る の で , 初期 の 要求 仕様 は 満た さ 
要求 精度 を 20 ビ ッ ト と し た うえ で , ブロ ッ ク 
演算 の けた の 割り 当て を 図 10 に 


に よる シミュ レ 
上 が 
な り ま す . 

図 を 図 9 に , 


, また RTL 記 述 を リス 


拡大 し て 示し ます . 


差 が 奇数 次 成分 に 


直線 補間 に 


よっ て 1 次 成 


分 を , 2 次 補正 に よっ て 2 次 成分 を それ ぞ れ 除 去 し た の で , 
3 次 成分 が 支配 的 に な る の は 先 の 級数 展 


を この 範囲 に 


見 立て る わけ で す . そし て , 図 7 を 正規 化 さ 


も な く 貞昌 
する な ら ば 


0.0000015 


0.0000010 


0.0000005 


較 
抽 


- 0.0000005 


- 0.0000010 時 


- 0.000001 


2 2.0 


8 2 次 補正 に よる 逆数 演算 の 誤差 
精度 が 20 ビッ ト 弱 ま で 上 が り , 初期 の 要求 仕様 は 満た され た こと に な る . 


開 の 例 を 出す まで 


理解 で きる で し ょ う . さら な る 精度 を 期待 
次 は 3 次 成分 を 除去 すれ ば よい の で す が , 演 


注 1: 


15: 11 
! 10: 0] 


! 15: 11] 


! 10: 0] 


和み" 0 で マク ロー リン 


展開 で き な い . 


し か し , 式 1) の よ 


に 定義 域 を 決め て いる の で , 本 1( 1 x) と し て 定義 域 を 0 ミ 
ェ < ご 1 と し て や れ ば , 級数 展開 も 可 台 
ン 展 開 な ど ( 


は 局所 的 に 


に な る . 


し か 精度 3 香 ら れ な い の で , 


始点 較 


JP 区 間 終 点 図 
テー ブル 区 


ゲイ ン ・ 図 
テー ブル 凶 
正規 化 
2 次 誤差 鐘 
補正 関数 図 


次 吉 差 
間 2 次 大 人 メ メン 


図 9 2 次 補正 追加 に よる 逆数 演算 ブロ ッ ク 


直線 補間 に よる 回 路 に 2 次 誤差 補正 機能 を 追加 し , 


16 12 8 4 0 
9 [1 0.000002 
i20 16 12 8 4 0 
{ag), ( a+ 万 層 
20 16 12 8 4 0 999909 
直線 補間 値 図 
上 12 ビ ッ ト 図 。 夫 
2 次 誤差 ゲイ ン 図 9 
正規 化 2 次 誤差 較 
ー コ - 0.000001 
2 次 補正 値 鐘 : 
20 16 12 8 4 0 
9 ー- 0.000002 
図 10 2 次 補正 追加 に よる 逆数 演算 の けた 割り 当て 図 11 


精度 向上 の た め , けた を 追加 し た . 


132 Design Wave Magazine 2006 Jonuary 


ーー つ 


も っ と も , マク ロー リ 
汎用 性 は な い . 


正人 図 


演算 精度 を 向上 させ る . 


9 


2 次 補正 に よる 最初 の 区 間 の 誤差 
区 間 内 の 誤差 が 奇数 次 成分 に よっ て 占め られ て いる こと が わか る . 


c C c 8 c 
C C e C に @ C C C c 
C C C ド C 

算 の 次 数 が 高く な る と 当然 , 演算 回 路 が か さば っ て きま す . 


な に か , 逃げ 道 は な い の で し ょ うか . 
@ さら に 2 次 以上 の 補正 を テー ブル 化 

葉 10) は 2 次 関数 な の で , これ で 補正 を か ける と 2 次 成 
分 に し か 効か な い の は 明らか で す . せっ か く テー ブル を 利 
用 し て いる の で すか ら , 2 次 以上 の 成分 を まとめ て テー ブ 
ル 化 し て , 式 10) の 代わ り に 使用 すれ ば , 高 次 成分 も 除去 
で きる は ず で す . この 方 法 を 検討 し て み ま し ょ う . 

ここ で ,「 だ っ た ら 1 次 成分 も テー ブル に 持ち 込ん で , ま 
と め て 補正 し て し まえ ば よい の で は 」 と 考え る 方 も いる か 
も し れ ま せん . 考え か た と し て は まち が っ て いな い の で す 
が , 1 次 成分 は 2 次 以上 の 成分 と 比べ 非常 に レベ ル が 高い の 


リス ト 2 図 9 の 回 路 の RTL 記述 


modu1e diy table b ( 
, // quotienE [20] : int [19:00] : Eraction 
d ): // dQdivigor [15 : 00] : Eraction ( [16] : statio in 1) 


// io 
npu モ も [15:00] d : 
outpu も [20:00] qd : 


// function (演算 処理 部 ) 
// tab1e 
Funot1on 


Ftab1e 
// [52:311] : 始 点 [30:09] : 終 点 [08 : 00] : ゲ イン 


[52:00] 


npu モ [15 : 11] d 
ご で a8@6 ( d ) 
5!h00 H Etab1e 


5!h02 8 Etab1e 
中略 
5 !h1d H Etab1e 
5!h1e Etab1e 
5'haf Etab1e 
endoage 
endFunot1ion 


0.00000015 


0.00000010 


0.00000005 


- 0.00000005 


- 0.00000010 


- 0.00000015 
l 旧 1.2 1.4 1.6 1.8 2.0 


図 12 高 次 補 補 正 に よる 逆数 演算 の 誤差 
23 ビ ッ ト 層 の 精度 を 得る こと が 確認 で きる が , 期待 され た 飛躍 的 な 精度 の 向 
上 に は 至っ て いな い . 


= 人 22'h200000,22「h1F07C2,9'!h1E9 } : wire [20:00 
5'h01  : Etab1e = { 22'h1F07C2,22'h1E1E1E,9'h1BE ) : 
= 人 22'h1E1E1E, 22「h1D41D4,9'h199 } : wire [17:00 


で , テー ブル 自体 が と て も 大 きく な っ て し まう の で すま 2. 
補正 テー ブル を 区 間 ご と に 持て ば , 非常 に 高い 精度 を 得 
られ ます が , 1 区 間 の アド レス だ け で も 11 ビ ッ ト ある の で , 
これ を すべ て の 区 間 に 持 つこ と は 避け る こと に し まし ょ う . 
1 区 間 の 正規 化 さ れ た 代表 値 を テー ブル 化し , 先 の 2 次 補 
正 の 場合 と 同じ よう に 区 間 ご と の ゲイ ン を 変え て , 直線 補 
間 値 か ら 減算 する こと に し ます . 代表 値 に は , 直線 補間 に 
よる 最初 の 区 間 の 誤差 を 正規 化し て 用 いま す . 

この 補正 方 式 に よる 誤差 を 図 12 Excel に よる シミ ュ レ 
ーション 結果 ) に , ブロ ッ ク 図 を 図 13 に , 演算 の けた の 割 
り 当て を 図 14 に 示し ます . な お , リス ト 1 と リス ト 2 に 基 
本 的 な 構造 が 示さ れ て いる の で , この 回 路 の RTL 記 述 は 
割愛 し ます . 

図 12 か ら , 28 ビ ッ ト 弱 の 精度 を 得る こと が 確認 で きま 


注 2: 2 次 まで は 演算 で 補正 を か け て 3 次 以上 を テー ブル で 処理 し 
た ほう が , テー ブル を 小さ く で きる . 本 稿 で は , 演算 回 路 を 
肥大 化 さ せな いと いう 論旨 で ここ に た どり 着い て いる ので, 
2 次 成分 は 含め る こと に する . 


Wire Pp_St : 

Wire p_end : 

Wire : gain : 

asg1gn {p_stt,p_end, gatn} 

= Etable (d[15:11] ) 


// interpo1ation 

Wi エ G 33:00 ntrpa = ( 12'h800 - d[10:00] ) 

* p gt 上 + d[10:00] * p end : // 直線 補間 
err2g 上 = ( 12'h800 - qd[10:00] ) 

* gd[10:00] + 11'h400 :  // 正規 化 2 次 関数 
1ntrpb = err2s 上 [20:11] * gann : 


// 2 次 誤差 生成 


Wire 21:00 1ntrD 


{ 22'h10C971,22'h108421,9'h048 } : = ntrpa [33 :11] - intrpb[17:09] + 1'b1 : 
{ 22!h108421,22'h104104,9'h045 H 
{ 22!h104104,22'h100000,9'h042 


// 直線 補間 - 2 次 誤差 


a881qdn fe = 1ntrp [21:01] : 


endmodu1e 


J 層 [ 15: 11] 図 

' 直線 図 

「15: 01 図 も 
高 次 較 
正人 


区 間 誤差 較 
ゲ ィ イン ・ 図 


図 13 高 次 補正 追加 に よる 逆数 演算 ブロ ッ ク 
2 次 関数 に よる 補正 回 路 の 代わ り に , テー ブル に よる 補正 回 路 を 追加 し , 精 
度 を 向上 させ て いる . 
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す が , 期待 され た 飛躍 的 な 精度 の 向上 に は 至り ませ ん で し 
た .「 すべ て の 次 数 を 含ん で いる テー ブル を 用 いた の だ か 
ら , 演算 誤差 以外 は 除 示 され る の で は な いか 」 と 期待 され 
た 方 も いる で し ょ う . この 原因 は , 直線 補間 に よる 最初 の 
区 間 の 誤差 を , 代表 値 と し て すべ て の 区 間 に 用 いた こと に 
あり ます . 実際 に は 各区 間 に お ける 高 次 成分 の 比 は 区 間 ご 
と に 違っ て いる の で す . し た が っ て , ゲイ ン を 区 間 ご と に 
調整 し て も , その ゲイ ン は 支配 的 な 2 次 成分 に 合わ せる こ 
と に な る た め , 3 次 以上 の 成分 に 対し て は 不 完全 な 補正 と 
な っ て し まう の で す . 最初 の 区 間 は 見 事 に 誤差 が 除去 で き 
て いる こと , その ほか の 区 間 で は 3 次 成分 が 支配 的 で ある 
こと か ら , その よう す が よ く わか り ま す . 


9 |1 


24 20 16 12 8 4 


{ag), 
{ ag+ 革 


直線 図 i24 20 


補間 値 


高 次 誤差 


ゲイ ン 


代表 


高 次 誤差 


高 次 


較 図 図 較 図 図 較 ) 


補正 値 
20 16 12 8 4 


IN 
選 - 


9 


14 高 次 補正 追加 に よる 逆数 演算 の けた 割り 当て 
精度 向上 の た め , けた を 追加 し た . 


テー ブル を 引く 基準 位置 凶 
fd 15: 111]) 図 


補間 位置 


id 10: 0] 図 
15 B_Spline に よる 補間 原理 


ディ ジタル ・ フ ィ ル タ と 同じ よう に 畳み 込み を 行う . 1 区 間 の 分 解 能 分 の 係 
数 が 必要 . 
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人 @ 補間 関数 の 高 次 化 
ここ まで は , 最初 に 直線 補間 あり き で , そこ で 取り きれ 
な か っ た 高 次 成分 を どう 除去 する か と いう 視点 で 考え て き 
まし た . 逆数 関数 は 曲線 な の で , 直線 補間 で は 取り きれ な 
い 成 分 が か な ら ず 付き まとっ て いた か ら で す . 次 は , 少し 
視点 を 変え て み ま し ょ う . 
ディ ジタル 信号 処理 に お ける 補間 処理 で は , 直線 補間 は 
きわ め て ベー シッ ク な 手段 で す . より 良い 結果 を 求め る 場 
合 は , 次 数 の 高い ディ ジタル ・ フ ィ ル タ が 使用 され ます . 
「 区 間 ご と の 離 散 的 な 代表 値 を それ らし く 補間 する 」 と いう 
視点 に 立て ば , 関数 の 波形 を 再現 する こと に も 応用 で き そ 
う で す . ここ で は , 補間 関数 と し て 3 次 の B_Spline を 採用 
し まし ょ う ( 右 掲 の コラ B _Spline 関数 」 を 参照 ). 
3 次 B_Spline の 補間 原理 を 図 15 に 示し ます . 図 15 に お 
いて , 補間 位置 で の 関数 値 は 式 11) で 表 さ れ ま す . 


3% 
7 の ⑦= う 2 ae 《 11) 


ん =0 凶 
この 式 11) の 演算 を 回 路 化す れ ば よい わけ で すか ら , こ 
こま で 使用 し た 区 間 の 始点 / 終 点 テ ー ブ ル の 代わ り に , g。 
て gs の 値 を 区間 ご と に テー ブル と し て 持つ こと に な り ま 
す . も っ と も と gs は , 区 間 の 始点 と 終点 と まっ た く 同 


区 間 分 割 に よる 図 
逆数 テー ブル 較 


[ 15: 11] 図 の 
d 凶 に | 9o テ ー ブ ル 較 
- 15: o]M 


[ 10: 0] 


す 


| 


1 区 間 内 図 
補間 係数 テー ブル 凶 


図 16 B_Spline 補間 に よる 逆数 演算 ブロ ッ ク 図 


9 の 上 位 で 区 間 ご と の 逆数 テー ブル 4 点 を 引き , 下位 で 係数 テー ブル を 引い 
て 畳み 込み を 行う . 


ーー 、 ーーー 生生 、 央 5 昼 回 計 計 8 
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ーー ーー C= ニ ニニ で = デー ーー ーー ーー es 且 
@ C C @ C "3 (7 し ーー /) ミー 
(ひめ J ーー だ 


じ も の で すし , go と gs は 1 区 間 隣 の 始点 と 終点 に な り ま す . 


どう いう テー ブル が 必要 か は あえ て 説明 する まで も な いで 990 
し ょ う . さら に , 新しい 要素 と し て B Spline 関数 の 発生 回 866me 
路 が 必要 に な り ま す . 数 値 演算 で 発生 させ る の も よい の で 
す が , 話 の 流れ か ら は テー ブル で の 表現 と な る で し ょ うか . oooooog 
この あたり は , 状況 に 応じ て 判断 し まし ょ う . ] 
16 に , ブロ ッ ク 図 を 示し ます . B_Spline 関数 は テー MM 作 M 衣 剛 
ブル で 表現 し て あり ます . 乗算 器 が 四 つ 並 ん で いる の が 回 川間 間 WW 
路 規模 の 面 か ら は 気 に な る と ころ で す が , 構造 が 単純 な の 1 時 re e9 
で 4 サイ クル か け て 1 系 統 で 済ま せる こと も 可能 で す . 誤 図 17 B_Spline 補間 に よる 逆数 演算 誤差 
差 の 評価 を 図 17 に 示し ます . 補正 な ど を 加え ず に , この 補正 な ど を 加え ず に , この 時 点 で 16 ビ ッ ト 弱 の 精度 が 確保 され て いる . 
時 点 で 16 ビ ッ ト 弱 の 精度 が 確保 され て いま す . また , 1 区 
間 内 の 誤差 に 着目 すれ ば , 1 次 成分 が 支配 的 な の は 明白 で また , これ ら の テー ブル 手法 が 逆数 関数 に 限ら れ た も の 
すか ら , 補助 的 な 補正 を か ける こと で 精度 の 向上 が 十分 期 で は な いこ と は すでに お 気づき で し ょ う . テー ブル を ROM 
待 で きま す . で な く RAM で 実現 し , 外部 か ら 内 容 を 書き 換え る こと で 
ネ ネ ネ ほか の 関数 も 実現 で きま ず 関数 の 特性 が それ ぞ れ 違う の 
今回 紹介 し た テー ブル に よる 手法 は , あく まで 1 例 と 考 で , 若干 冗長 性 が 必要 に な る と 思う が …). 
えて くだ さい . 要求 され る 精度 , リソー ス , 遅延 レイ テ 
ン シ ) な ど 状 況 に 応じ て , 実現 する 形 は どの よう に で も な 
る と いう こと で ず テー ブル と 演算 の 分 担 な ど は , チュ ー 9 すさ Le っ 
、 スタ ンダ ー ド ・ ロ ジッ パ 株 ) 
ニン グ に よっ て 回 路 規模 に 大 き な 影 響 を 与え る る). また , ほ 
か の 手法 を 組み 合わ せる こと で お も し ろ い 結果 が 得 ら れる べく 筆 者 プロ フィ ー ル > 
、 鈴木 昌治 . 枯れ た 回 路 と 思わ れる 部 分 で も , 何 か お も し ろく な ら な い 
可能 性 も あり ます . 柔軟 な 発想 で 新しい 実現 方 法 を 模索 し 才 0 上 還 史 挟 12 の eS まで の の 本村 
て みる の も 楽し いで し ょ う . に ら ま れ な い 範 囲 で 」 で す が … 


B_Spline 関数 


Splin& スプ ライ ン ) 関数 は , 離散 的 な 点 を 適当 な 滑ら ちか さ で 連結 図 B-1 は , その 関数 波形 を 示し て いま す . 
させ る 関数 で す . 平面 上 の 離散 点 を 自在 定規 で 結ぶ 操作 を 数 学 的 に 式 B-1) が 与え られ て いれ ば , 補間 係数 テー ブル を 作成 で きる の で , 
解決 する 方 法 と 考え る こと が で きま す . スプ ライ ン 関 数 に つい て 詳細 を 知っ て いな く て も 関数 の 近似 を 行う こ 

スプ ライ ン 関 数 の 用 途 は 多岐 に 渡り ます が , 関数 の 近似 も 有効 な と が で きま す . も ちろ ん , テー ブル を 用 い ず に 演算 で 補間 係数 を 求 
アプ リケーション の 一 つ で す . 本 稿 で は , テー ブル 上 の 各区 間 の 境 め る こと も で きま すね . 
界 点 を それ らし く 結ぶ こと で 区 間 内 を 補間 し , 関数 を 近似 する た め 
に 使用 し て いま す . 

B_Spline は , スプ ライ ン 関 数 の 中 で も っ と も 基本 的 basic) な 関数 
で あり , 3 次 の 場合 は 以下 の よう に 定義 され ます . 
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図 B-1 
B_Spline 関数 
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